Ground States of the SU(A^) Heisenberg Model 
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The SU(A'') Heisenberg model with various single-row representations is investigated by quantum 
Monte Carlo simulations. While the zero-temperature phase boundary agrees qualitatively with the 
theoretical predictions based on the expansion, some unexpected features are also observed. 
For N > 5 with the fundamental representation, for example, it is suggested that the ground states 
possess exact or approximate U(l) degeneracy. In addition, for the representation of Young tableau 
with more than one column, the ground state shows no VBS order even at A*' greater than the 
threshold value. 

PACS numbers: 75.40.Mg, 75.10.Jm 
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Since the resonating-valence-bond (RVB) state was 
proposed[l| as a possible mechanism that supports novel 
super conductivity in cupurates, it has been a major tar- 
get of condensed matter theory to find a short-range in- 
teraction model that realizes a spin liquid state at zero 
temperature. Introducing frustrations into the model 
that otherwise has a magnetic ground state is one of 
promising directions to achieve this goal. Numerical ver- 
ification for such models is, however, technically very 
difhcult, and conclusive results are still missing for ge- 
ometrically frustrated systems such as the antiferromag- 
net on a triangular lattice. Another approach was taken 
in [1, where the authors generalized the Heisenberg 
antiferromagnet to higher symmetries, thereby increas- 
ing the model's degrees of freedom and enhancing quan- 
tum fluctuations. Based on the 1/7V expansion treat- 
ment, they predicted that the model with sufhciently 
large A'' has a valence-bond-solid (VBS) ground state 
with spontaneously-broken lattice symmetry. The nature 
of the ground states may depend on the representation, 
somewhat analogous to what happens in the SU(2) mod- 
els in one dimension though the underlying mechanisms 
may be rather different. 

More specifically, it was sugg estedH H that, for the 
model with the Young tableaux with m rows and n 
columns, the N-n phase diagram does not strongly de- 
pend on m, and has a line of phase transition separating 
the small- Neel region from the large- VBS region. It 
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FIG._b Three types of VBS states: (a ) columnar {D^ / 
and Dy — 0), (b) nematic {D^ = Dy — 0), (c) plaquette 
(D^ — Dy 7^ 0). The quantities and Dy characterize 
spontaneous translational-symmetry-breaking in x and y di- 
rections, respectively. (See Eq.Q for the definition.) 



was also argued that the nature of the VBS ground state 
can be classified according to the quotient of the division 
of 71 by 4. If n = 1 or 3 (mod 4), the ground state has 
a columnar ordering (Fig.[T] (a)) with the translational 
symmetry and the 90-degree rotational symmetry both 
broken, whereas if n = 2 (mod 4) it has a "nematic" 
VBS ordering (Fig.[T](b)) with only the lattice-rotational 
symmetry broken. Finally, if n is a multiple of 4, there 
is no spontaneous breaking of the lattice symmetry. 

They also showed that the 0{1/N) effective Hamilto- 
nian with the fundamental representation is exactly the 
same as the quantum dimer model |^ at = 0. For 
the quantum dimer model, it was shown by a numerical 
calculation^ that the ground state is a VBS state with 
columnar arrangement of dimers. This is in agreement 
with the above conjecture that the degeneracy is 4 for 
n = 1 . A direct check of the spontaneous breaking-down 
of the translational symmetry for m = n = 1 was carried 
out in the previous paper , which yielded the transition 
value of A^, namely 4 < iV*(m = n = 1) < 5. 

The model Hamiltonian we discuss is defined as 
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H = j J2 



(1) 
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where the operator S"^ is the generator of the SU(A^) 
algebra. Here we consider the simple square lattice with 
the periodic boundary condition. We divide the whole 
lattice into to two sublattices, say, A and B. The rep- 
resentation of the generators on sublattice A is charac- 
terized by the Young tableau with a single raw (m = 1) 
and varying number of columns. The representation on 
sublattice B is the conjugate of that on sublattice A. 
We have performed quantum Monte Carlo simulations 
by the directed loop algorithm|71] and extrapolated the 
results to the zero temperature limit. The algorithm can 
be obtained by generalizing the idea of the coarse-grained 
algorithm for A^ = 2 to the present general A^ case. 
The technical details will be published elsewhere. Q 

Our principal findings are as follows: (i) For n = 1, 
the ground state in the large-A^ region adjacent to the 
Neel region is the columnar VBS state as the previous 
works suggested. However, up to L = 32, it appears to 



2 





L=4 B 




L=8 ■ 




L=16 e - 




• L=32 • 




A L=64 6 







2 3 4 5 

N 



U 



10" 

10-^ 
10-2 
10"^ 
10-^ 
10^ 





1 

N= 4 e ■ 




N= 5 A 




N= 6 ^ _ 




N=15 K : 


1 1 1 


I 

1 



16 



32 



FIG. 2: The correlation ratio, Rs{L), of the magnetization 
for the fundamental representation {n — 1). 

have infinite degeneracy corresponding to a continuous 
U(l) symmetry that does not exist in the original micro- 
scopic Hamiltonian but may emerge asymptotically. Ac- 
cordingly, the resonons [3| are gapless. (ii) The ground 
state in the large- iV region for n > 1 shows no evidence 
for break down of any lattice translational symmetry, in 
particular, the ones predicted in Read and SachdevjH, 
The distance dependence of the energy-energy correlation 
function is consistent with the algebraic decay. 

In the present letter, we start from showing some re- 
sults that confirm the conclusions of the previous nu- 
merical work on the SU(A^) model, and then discuss the 
new findings. We first look at the case of the funda- 
mental representation {n = 1). It was found 0] that 
for TV < 4 the ground state is a Neel state whereas it 
does not have a spontaneous staggered magnetization for 
N > 5. Instead the ground state for = 5 and 6 was 
found to possess the spontaneous dimerization. Although 
a macroscopic quantity is often used for a probe for de- 
tecting the transition, here we examine two-point correla- 
tion functions in order to see the transition more clearly, 
Cq{R;L) = (g(i?e,0Q(O)}L-(g(i?e^))(g(O))L where i 
is the system size, (/i = x, ?/) a lattice unit vector, and 
Q{R) an arbitrary quantity locally defined around the 
position R. We also use the correlation ratio [l^] defined 
simply as Rq{L) = Cq{L/2; L)/Cq{L/4:; L). Similar to 
the Binder parameter, the correlation ratio is a dimen- 
sionless quantity and, when plotted against a relevant 
physical parameter, the common crossing point of curves 
with various system sizes marks the transition point. In 
Fig. [21 the correlation ratio for the "magnetic" moment 
is plotted against A^ for various system sizes for to = 1. 
Here, the "magnetic" moment is defined as 

S(fi) = 5"(H)- 522(H). 

While the ratio tends to converge to unity for A^ = 2, 3 
and 4, indicating the existence of the Neel ordering, it 
rapidly decreases for N > 5. Based on this figure, we 
can estimate the transition point as N*{n = 1) ~ 4.3. 

While Fig. [5] establishes the absence of the magnetic 
ordering for A^ = 5,6,7, it does not tell us much about 



FIG. 3: The two-point correlation function the rotational- 
symmetry breaking order parameter. A, for the model with 
the fundamental representation {n — 1). 

the nature of the ground state. It was shown in that 
in the ground state the lattice translational invariance 
is broken. Here we reconfirm the presence of the VBS 
order by detecting the long-range correlation in Cyi(iJ; L) 
where A is defined as 

A{R)=P^{R)~Py{R), (2) 

and P^{R) is the nearest-neighbor product of the "mag- 
netic" moments, 

Pf,{R) = j:{R)^{R + e^). (3) 

In Fig. [21 Ca{L/2; L) is plotted. It is expected that the 
correlation must decay down to zero in the large L limit 
for a Neel state, whereas for any type of the VBS states 
in Fig.[l] it must converge to a non-zero value. Fig. [31 
clearly shows that the large- A^ phase for the fundamental 
representation is the VBS phase. 

These observations are consistent with the columnar 
VBS states (such as (a) in Fig.[T]) as is conjectured for 
N > N* [1, y]- However, the above-mentioned probe A 
only is not sufficient to rule out the plaquette VBS state 
((c) in Fig.[T]). Since the columnar VBS states do not 
possess the 90-degree rotational symmetry while the pla- 
quette VBS state does, one may naively expect that mea- 
suring the difference between the average bond strength 
in the x direction and that in the y direction can distin- 
guish the two. Therefore, we have computed the quan- 
tity B = {[Px — Py) ) where an overline indicates the 
average over the volume, i.e., P^i = '^rPi^X^)- The 
quantity B is plotted in Fig. [31 against the system size. 

Instead of distinguishing the two types of states, the 
quantity B turns out to reveal an interesting property. 
The quantity B is proportional to l/L^ not only for 
A^ = 3 and 4 but also for N = 5 and 6, indicating that 
the average bond strength in the x direction is the same 
as that in the y direction even in the VBS state. It fol- 
lows that the expectation value of the nearest-neighbor 
correlation P^i{R) in a single (pure) VBS ground state 
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FIG. 4: The squared difference of the average bond strengths 
in the x and y direction for the fundamental representation. 
The straight hues that correspond to B oc are drawn for 
comparison. 

has the form 

(P^(i2))si„gie = P + :D;x (-1)^- {ii = x,y) (4) 

with the direction-independent average value P and 
some direction-dependent constant that char- 

acterizes the dimerization order. The constants 
Dx and Dy are expressed by the local quantities 
D^{R) &sW^^ V-^Y.RDf.{R) where D^{R) = 
{P^{R + e^)-P^{R))l2. 

In order to distinguish the columnar state from the 
plaquette state, we have to examine the joint distribution 
function, Prob(_D2;, Dy). For small systems, it consists of 
a single broad peak at the center, {Dx,Dy) — (0, 0), even 
for TV > 5. However, as the system grows the weight at 
center should diminish and peaks must appear according 
to the structure of the VBS phase. If the columnar VBS 
states are the true ground state, four peaks must develop 
at symmetric positions on x and y axes ( {Dx,Dy) — 
(±Z?, 0), (0, ±_D) ) whereas if the plaquette VBS states 
are the ground states the peaks must appear at diagonal 
positions ( (D^,D^) = ±{D,D),±{D,-D) ). 

As we see in Fig. [SI the distribution obtained from our 
computation shows neither feature. The distribution is 
circularly symmetric. This feature does not depend on 
N, the temperature, or the system size (at least up to 
L ~ 32), whenever a finite VBS order is observed. This 
suggests that the ground state is not only four-fold degen- 
erate but infinitely degenerate. The ground state mani- 
fold possesses the U(l) invariance although it is not obvi- 
ous from the microscopic Hamiltonian Eq. ([T|) . Since each 
ground state breaks this U(l) symmetry, there must be 
a Goldstone mode associated with it, which may corre- 
spond to "resonons" mentioned in Q, though its gapless 
nature was predicted only for an isolated critical point. 
If the apparent U(l) symmetry persists to the thermody- 
namic limit, the resonons must be gapless not only at the 
isolated critical point N — N* but also in the whole re- 
gion of the VBS phase. However, another possibility can- 
not be excluded by the present numerical calculation, i.e.. 
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FIG. 5: The frequency of observing a pair of values {Dx, Dy) 
during the Monte-Carlo simulation for A*' = 6 at L = 32 and 
l3 = 16. The value at each pixel is the average over eight sym- 
metrically equivalent points, {±Dx, iDy) and {±Dy,itDx) 
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FIG. 6: The correlation ratio of the magnetization for the 
representation of two-box Young tableau (n = 2). 



the possibility that this apparent U(l) symmetry may be 
a transient behavior that applies only to certain inter- 
mediate length scales and the discrete symmetry of the 
original microscopic Hamiltonian is recovered for larger 
length scales. If this is the case, there should be a fi- 
nite number of preferred angles in the two-component 
order-parameter space to which the true ground states 
correspond, and the resonon modes have a small but fi- 
nite excitation gap. 

Next we consider the model with the representation 
of Young's tableaux with two or more columns {n > 2). 
In order to determine the critical value of N, we have 
computed the correlation ratio for the magnetization as 
described above. For n — 2, the correlation ratio for 
the magnetization is plotted in Fig.[Sl The Neel state 
is the ground state for N < 9 whereas it is not for 
N = 10 and larger. We can see the trend changes 
around N*{n ~ 2) ^ 9.5 at which curves cross each 
other. This is again in a good agreement with 1/A'^ 
expansion result, N* ~ 5.3 n 0, For n = 3 and 

n = 4 we have performed similar computation and lo- 
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FIG. 7: The zero-temperature phase diagram of SU(A'') model 
on a square lattice with single-row (m = 1) representations. 
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FIG. 8: The correlation function of the rotational-symmetry- 
breaking order-parameter A at the largest distance for n — 
1, 2, 3 and 4. The value of A'' for each n is close to but larger 
than the ythreshold value. 



cated the the phase boundary as N*{n = 3) 14 and 
N*{n = 4) ^ 20, respectively. These results are summa- 
rized in the schematic phase diagram (Fig. [7]). 

Lastly, we turn our attention to the nature of the non- 
magnetic ground states for n > 2. As mentioned above, 
any type of the VBS states in Fig. [T] can be characterized 
by a non-zero value of Ca{L/2; L) in the L — > oo limit. 
This correlation function is plotted in Fig.[8]for n = 1,2,3 
and 4. Characteristic features do not depend on N for 
each n. We only show the data for a single value of N 
for the sake of readability of the figure. (The value of 



N shown in Fig. [8] is chosen so that it is close to but 
definitely above the estimated threshold value.) The top 
left panel for n = 1, iV = 5 is shown for comparison. As 
is clear from the figure, the system possesses the VBS 
ordering only for n — 1, but not for n — 2, 3, 4. 

To obtain some information as to whether the large- iV 
phase at n > 2 is gapless or not, we have also com- 
puted the correlation function of P^{R) itself, which 
is proportional to the energy-energy correlation. In a 
clear contrast to the case of n = 1, we have seen mono- 
tonic decrease to zero for n — 2,3,4 which is consistent 
with the absence of the rotational symmetry breaking for 
n = 2,3,4 mentioned above. We do not clearly observe 
exponential decay in any of the disordered cases. The 
decay in the correlation up to the length scale of L ~ 32 
is consistent with the power law with the decay exponent 
between 2 and 3. However, since the correlation function 
for large L is very small and the relative error is too large 
to obtain a definite answer to the question whether this 
phase is gapless or not. 

To summarize, in addition to the confirmation of the 
previous result ,6l|, we have found an evidence for an 
emergent (exact or approximate) U(l) symmetry of the 
ground state space of the SU(N) model with the funda- 
mental representation, and also found no VBS order in 
the large N region adjacent to the phase boundary to 
the Neel region. The latter finding is the first strong evi- 
dence for the spin-liquid (disordered) ground state in the 
present model that does not have geometrical frustrations 
in contrast to the models studied previously in search 
for the quantum disordered states. As for the first find- 
ing, it is appropriate to mention the recent works on the 
deconfinement critical phenomena (DCP). The present 
SU(jV) rnodel is the model that was discussed in previous 
works[l3, 0, 13 related to DCP, where the authors ar- 
gured that some unconventional type of the second-order 
phase transition is possible between two phases with ap- 
parently unrelated symmetries (e.g., the VBS phase and 
the Neel phase). There also discussed was the emergent 
U(l) symmetry at (or near[l3|) the critical point. We 
suspect that the U(l) structure that we have observed 
reflects validity or approximate validity of the DCP sce- 
nario in the present model at zero temperature. 

The author would like to thank Cristian Batista and 
Anders Sandvik for useful comments. The computation 
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Supercomputer Center, Institute of Solid State Physics, 
University of Tokyo. 
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